Comment on "Do Earthquakes Exhibit Self-Organized Criticality?" 

In a recent Letter, Yang, Du, and Ma 111 study the interesting problem of the temporal structure 
of seismicity and its relation with self-organized criticality (SOC). Their main finding is that the 
reshuffling of earthquake magnitudes changes the shape of the earthquake recurrence-time (or 
first-return-time) distribution when the low-magnitude bound, M c , is raised. Subsequently, they 
conclude that it is not true that an earthquake cannot "know" how large it will become. First, we 
show that this important implication is unjustified. 

Yang et al. have in mind a fully uncorrelated temporal point process with independent magni- 
tudes as a picture of SOC systems. It is obvious, by construction, that this model is invariant under 
random rearrangements of the data; as Yang et al. do not find this invariance in Southern California 
seisimicity they claim that "earthquakes do not happen with completely random magnitudes" and 
therefore they are not a SOC phenomenon. In fact, the only conclusion that can be drawn from this 
is that the seismicity time series is not uncorrelated, and there exists some dependence between 
magnitudes and recurrence times. [This conclusion can be obtained directly, from the fact that a 
scaling law exists for the recurrence-time distributions corresponding to different low-magnitude 
bounds, with a scaling function that is not a decreasing exponential ||2|] (characteristic of a Poisson 
process, the only uncorrelated process which verifies a scaling law).] 

The existence of correlations means that, for a given event /, its magnitude M; may depend 
on the magnitude of the previous event, Af,-_i, as well as on the backwards recurrence time, 7] = 
U — U-i, with ti and t{-\ the time of occurrence of both events. This dependence can be extended 
to previous magnitudes and recurrence times, 7J_i,Mj_2,7/-2> etc. But further, the recurrence 
time to the next event, 7}+i, may depend on the previous magnitudes, Mj and recurrence times 
Tj, j < i. The reshuffling of magnitudes performed in Ref. yj] breaks (if they exist) the possible 
correlations of M ; with the previous magnitudes, as well as with the previous recurrence times, and 
the correlations of T i+ i with the previous magnitudes (but not with the previous recurrence times). 
Therefore, any of the influences M{-\ — > Mj, 7} — > Mj, or Mj — ► 7]+i, may be responsible of Yang 
et a/.'s results. 

The most direct way to test the dependence of a given variable, in this case Mj, with an- 
other variable X, is to measure the probability density of X conditioned to different values of 
Mi, P(X\Mi), and compare with the unconditioned probability density of X, P(X). This is what 
Fig. 1(a) displays, using X = 7j and X = T i+ \ for the same data as Ref. [1], but restricted to pe- 
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riods of stationary seismicity (otherwise, for strong aftershock sequences the recurrence times are 
shorter and more sensitive to catalog incompleteness). As P{Ti\Mj) remains practically unchanged 
for different sets of values of M ; , temporal causality leads to the conclusion that A/,- is independent 
on T(. In contrast, 7} + i clearly depends on M ; , as P(Ti + \\Mi) changes for different sets of values 
of Mj. In other words, the larger the magnitude Mi, the shorter the time to the next event 7}+i, but 
the value of this time has no influence on the magnitude of the event, M,-+i. On the other hand, 
Fig. 1(b) shows that P(Mi\Mi-\) turns out to be indistinguishable from P(Mj), ensuring the inde- 
pendence of Mi and M,-_i, V/ if the 7/'s are restricted to be larger than 33 min (shorter periods of 
time are not reliable, due to data incompleteness). So, when an earthquake starts, its magnitude 
is undetermined (at least from the information available at the catalogs), whereas the time to the 
next event decreases when that magnitude turns out to be large. 

A second, independent point to clarify is the identification of SOC with the total absence of 
correlations. It is true that the BTW sandpile model displays an exponential distribution of recur- 
rence times, but SOC is much more diverse than the BTW model. For instance, the Bak-Sneppen 
model or the Oslo-ricepile model are two well recognized examples of SOC with totally different 
recurrence-time distributions. Finally, it is necessary to stress that the concept of SOC (as it hap- 
pens with chaos) does not exclude the possibility of prediction, as Ref. [17] of Yang et al. clearly 
showed. So, nothing in Ref. yj] is against the SOC picture of earthquakes. 
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FIG. 1: (color online) (a) Probability densities P(7}|M,-_i) and P{Tj\Mi) (shifted upwards) compared to 
P(Ti). (b) Probability density P(M,-|M;_i) compared to P(M,) with T t > 2000 s. 
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